Decomposition of force fluctuations far from equilibrium 
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By studying a nonequilibrium Langevin system, we find that a simple condition determines the 
decomposition of the coarse-grained force into a dissipative force, an effective driving force and noise. 
From this condition, we derive a new universal inequality, D > 7/idT, relating the diffusion constant 
D, the differential mobility /xa, the bare friction constant 7 and the temperature T. Due to the 
general nature of the argument we present, we believe that our idea concerning this decomposition 
can be applied to a wide class of systems far from equilibrium. 
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The nature of a force depends on the scale on which it 
is observed. For example, consider a force exerted by wa- 
ter molecules on a colloidal particle. Such a force can be 
described by mechanical laws on molecular time scales, 
while it is described as a dissipative (frictional) force and 
thermal noise on time scales of the order of 10 -3 sec. 
In analogy to this example, for a wide range of systems, 
including bio-mechanical systems and granular sys- 
tems 0, it might be expected that a fluctuating force 
obtained through some coarse-graining procedure can be 
decomposed into a dissipative force and other compo- 
nents. In the regime near equilibrium, such a decompo- 
sition is uniquely determined by fluctuation-dissipation 
relations (FDRs). However, no rule is known that deter- 
mines such a decomposition far from equilibrium Q . We 
wish to discover a rule of this kind for a class of systems 
exhibiting fluctuating forces. 

With the above-stated purpose, in the present work, 
we study a Langevin equation describing the motion of 
a Brownian particle with a tilted periodic potential in a 
one-dimensional space. Although we study the simplest 
system realizing nonequilibrium steady states (NESSs), 
the arguments below can be applied to a wide class of 
Brownian motors 0, |H E| . The Langevin equation that 
we analyze is described by 



(1) 



Here, 7 is a friction constant, U(x) is a periodic potential 
of period £, f is a constant external driving force, and £ 
is Gaussian white noise satisfying 



(t(t)Z(t'))=2 1 TS(t-t), 



(2) 



where T is the temperature of the environment, and the 
Boltzmann constant is set to unity. We consider a de- 
scription of large-scale motion, which is obtained by tak- 
ing an average over a time interval St that is chosen to 
be sufficiently longer than the characteristic time of the 
system. In this description, the finite time average of the 
force —dU/dx acts on the particle as a fluctuating force. 
We conjecture that this fluctuating force can be decom- 
posed into a dissipative force, an extra driving force and 
random noise. 



In order to investigate time averaged quantities, in- 
cluding that of —dU/dx, we introduce the finite time av- 
erage of an arbitrary quantity Z(t): 



— _ 1 

Zn = Tt 



dtZ(t), 



(3) 



tn 



where t n = nSt, n — 0, 1,2- •■. Then, the finite time 
average of —dU/dx that we consider is given by 



dU 

d,X n 



1 f dt 

St .L dx 



(4) 



C (t) 



We hypothesize that this can be decomposed into a dis- 
sipative component A(x n +i — x n )/St, where x n = x{t n ), 
and a non-dissipative component. That is, we assume 
the form 



dU j Xn+l - 

dx n St 



B n , 



(5) 



where A is a constant and B n is a fluctuating quantity 
whose statistical average (B n ) takes a non-zero value / p : 



Bn = f p + SB n 



(6) 



The quantities f p and SB n correspond to an extra driving 
force and random noise, respectively. Substituting ifKjl 
into an integrated form of we obtain 



(7 -A) 



St 



f + F + 6B n + Z n . 



Then, with the definitions 

r 

F 



j -A, 
f + P,_ 



we express {7} as 



r 



si 



= F + E n , 



(7) 



(8) 
(9) 

(10) 



(11) 



where H n is expected to exhibit a Gaussian distribution 
for large St. This equation is regarded as an effective 
model of QJ. 
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The main claim of the present Letter is that the simple 
condition 



st 



lim St (SB„i n ) = 



(12) 



uniquely determines the constant A in JSJ. Note that a 
correlation of time averaged quantities is proportional to 
St -1 in general and i|12[l indicates that the proportional 
constant for the case {SB n £ n ) becomes zero. The condi- 
tion l|12fl implies that A(x n +i —x n )/St, which fluctuates 
in time, can be distinguished from B n by the condition 
(|12fl lj. After presenting the proof of this claim, we re- 
mark on three important topics related to it: a general 
inequality obtained as a direct application of l|12|). ener- 
getic considerations related to (|12fl . and the relation be- 
tween (|12|) and a time reversal symmetry in the stochastic 
sense. 

Preliminary consideration: Before presenting the 
proof of the main claim, we first consider how the pa- 
rameters of the effective model (|ll|l can be expressed 
in terms of the steady state velocity v s and the diffusion 
constant D, defined as v s = lim^oo {(x(t) — x(0))/t) and 
D = lim^oo ((x(t) — x(0) — v s t) 2 /2t). Because v s and 
D are independent of the scale on which we describe the 
system, the same values of v s and D should be obtained 
from (fill) . This implies the relations F — Tv s and 



(S„S m ) = 2S nm DT 2 (StY 



(13) 



From these expressions, all the parameters of the effective 
model (lllf) are given in terms of v s and D when T is 
determined. 

In the equilibrium case (/ = 0), T should satisfy the 
relation DT 2 = FT, which is referred to as the FDR of 
the second kind [8j. From this, T is expressed as 



T 
D 



Furthermore, we can prove 



D = /zT, 

where [i is the mobility, defined as fi 
The two FDRs and JT5J lead to 



(14) 



(15) 



lim/^o v s (/)//. 



(16) 



However, for NESSs far from equilibrium, i|15|) is vio- 
lated ||, and l|14(l does not hold in general. Therefore, 
for treatment of such systems, it is necessary to find a 
guiding principle to determine T. 

In a previous work [9|, we studied NESSs and pro- 
posed a natural method to determine T by considering 
the response of the particle to a slowly varying potential 
V(x) in space. Below, we present a heuristic argument 
from which the result obtained there can be understood. 
For detailed presentation of the systematic perturbation 
method used to derive this result, see Ref. 

Because the gradient of the slowly varying potential 
V(x) can be regarded as a modulation of the external 



force /, the large-scale motion in the modulated system 
is described by 



X n +1 ~X„ I dV(x n ) 



St 



l — "n 



where /id is the differential mobility defined by 

dv s (f) 



df 



(18) 



Then, assuming that —dV/dx n is the force acting on the 
particle even in this effective description, from l|17l) we 
obtain the result 



(19) 



Note that this expression represents an extension of Ijlfijl 
to the presently considered non-equilibrium case. From 
© and (|19|) . we find that the constant A appearing in 
(JSJ should satisfy 



7 - A = M d 



(20) 



Proof: We demonstrate that the decomposition con- 
dition l|12fl uniquely determines the constant A, yielding 
(|2U|I . Using JSJ and Q), we rewrite the condition (|12H as 



i im 5t /( d E +A x ^\- Xr - 

St^co \ \ dx n St 



x n+ i - x n dU _ ( | 
St dx • 



(21) 



where ( } c represents the cumulant. Through the defini- 
tion 



G = - lim St(^ 

St— too \ dx n Ot 



it is easy to obtain 



7 - A = 



2 7 T 



(22) 



(23) 



2jD-G 

In order to connect lll'-'il) with the differential mobility 
m defined by l(TH)l . we express it in terms of D and v B . 
We start with the path integral representation 

Wit) " 1(0)1 -fw>-m) 

e-^rf-LH^-f+^T . (24) 

Then, differentiating both sides with respect to /, we 
derive 

| (x(St) - x(0)) = pt ^Lst + 0(St 2 ) (25) 

for large St. This leads to the relation 

2 7 D-G 



A*d = 



2 7 T 



(26) 



Comparing l|23|) and (|26|l . we have arrived at H2U|) . 
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Inequalities: As a simple application of the decompo- 
sition condition (|12|) . we derive several useful inequalities. 
From the square of both sides of pO(l and the condition 
(|12[) . we obtain 



2T 2 D = 6t ((6B n ) 2 ) + 2 7 T, 



(27) 



where we have used <|13[) . This immediately leads to an 
inequality relating the intensity of the force noise in the 
original system and the quantity representing its effective 
value in the coarse-grained system: 

T 2 D > 7 T. (28) 

Because V = /i^ 1 (see H190 ). I|28ll can be written as 

D > 7/ijT. (29) 

We believe that this inequality holds in other Brownian 
motors [H 0- E3] because the path integral expression is 
valid even for cases with a time dependent potential and 
the expressions <|21ll - (|26[) given in the proof are the same 
for those models. The inequality l|29(l involves only di- 
rectly measurable quantities and therefore can be tested 
exp er iment ally. 

In a related work, Sasaki conjectured that the inequal- 
ity D > /itjT holds generally for Brownian motors |lCj . 
If we define the effective temperature T e g using the FDR 
violation factor [9j, this conjecture is equivalent to the 
assertion that T e ff is not less than the temperature of 
the environment, T. Because it has been observed that 
T c g > T in glassy systems and driven many-body 
systems 0,^3), the inequality D > fi^T does seem plau- 
sible. However, Sasaki reported that this inequality is 
violated for the model with an appropriate choice of 
This result leads us to believe that, in the 
present context, if there exists a generally valid inequal- 
ity among measurable quantities, perhaps it involves the 
intensity of the force noise and its effective one, not the 
temperature and its effective one. 

Energetics: In equilibrium systems, heat is distin- 
guished from work according to the second law of ther- 
modynamics. However, obviously, heat can be considered 
as a mechanical work done by a force at a microscopic 
scale. This tempts us to investigate how heat and work 
come to be expressed in different ways through a coarse- 
graining procedure. We treat this problem on the basis 
of the decomposition of the force —dU/dx n given by JSJ. 

In Langevin systems, the heat absorbed from a heat 
bath is interpreted as the work done by a force —~/dx/dt+ 
£ exerted by the heat bath 0] . With this interpretation, 
the heat absorbed during an interval t n < t < i„+i can 
be expressed as 



dx 
'~dt 



n odx(t), 



(30) 



where the symbol o indicates that the integral here is 
the stochastic Stieltjes integral in the Stratonovich sense 



[l5j . Then, the energy balance equation for the model 
JTJ is derived as 

U(x n+ i) - U(x n ) = q n + f(x n+ i - x n ). (31) 

Through similar considerations applied to the effective 
model (|llf) , we define the heat (absorbed from an "effec- 
tive heat bath") as 

Qn = - {^ Xn+1 fi t ~ ( X n+1 - Xn)- (32) 

With this, the energy balance equation for l|ll(l is ob- 
tained as 



Q n + F(x n +i 



0. 



(33) 



The difference between the two quantities q n and Q n be- 
comes obvious when their steady state averages are com- 
pared; we have (q n ) = —fv s St and (Q n ) = —Fv s 5t, while 
it is known that outside the linear response regime, in 
general F ^ / (see Fig. 2 of Ref. @). 

In order to understand the difference between q n and 
Q n , we consider a decomposition of the work done by 
the force —dU/dx during a time interval t n < t < i„+i, 
which is written as 



dll dx 
-— — — of 

dx dt n 



(34) 



Although we conjecture that QP and W% correspond to 
"heat" and "work" , respectively, no rule is known that 
distinguishes "heat" from "work" in this case. However, 
in the present system, it seems natural to assume 



W%=P(x n+1 -x n ), 



(35) 



because the extra driving force / p was determined from 
the decomposition condition l|12|l. With this assumption, 
we can derive the relation 



Qn 



Qn' 



(36) 



which provides a clear interpretation of the difference be- 
tween q n and Q n . 

The argument above leads to the following question: Is 
there a simple rule of energetics from which we can ob- 
tain the decomposition l)34|) along with 135fl without the 
decomposition condition (|12|l ? This fundamental ques- 
tion will be studied in the future. In addition to their 
role in such fundamental problems, energetic considera- 
tions applied to different time scales may be useful when 
we attempt to interpret the efficiency of motor proteins 
[l(| within stochastic models [l7| . 

Adjoint dynamics: As a final topic here, we con- 
sider the relation between the decomposition condition 
(|12|l and time reversal symmetry. In order to rep- 
resent this symmetry explicitly, we consider the path 
probability density P for a discrete time series [x]n = 
[xq,x\, ■ ■ ■ ,xn) generated within the model under con- 
sideration. The time-reversed path probability density 
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P* is defined by P*([x] N ) = P([x] N ), where [x] N rep- 
resents the time reversed trajectory of [a;]jv, that is, 
x n = xn — n . When P* is obtained from the frequency 
distribution of trajectories in a steady state for some 
stochastic dynamics, these dynamics are called "adjoint 
dynamics" . 

When the constant A appearing in (JSJ is chosen cor- 
rectly, J7J can be regarded as an effective Langevin 
model. In this case, we find that adjoint dynamics are 
described by 

( 7 _ A fn + l^Xn = _f_ fP + SBn + (37) 

Then, note that the condition (|12fl is equivalent to 

(H5B„+0 2 )H(55 n +0 2 ). ( 38 ) 

Because both 8B n and £ n exhibit Gaussian distributions 
for sufficiently large St, the condition (|38|) allows us to 
replace 5B n + £ n in (p?7jl by — 5B n +£ n - Thus, the adjoint 
dynamics can be expressed as 

{l -A f n+1 - Xn =-f-B n +l n . (39) 

Then, because Q can be rewritten as 

in-A) Xn+ \~ Xn =f + B n +l n , (40) 



we find that the decomposition of —dU/dx n given by © 
is characterized by parity with respect to time reversal 
in the stochastic sense. More specifically, the dissipative 
force A(x n+ i —x n )/St remains, while the other contribu- 
tion, B n , changes sign in the adjoint dynamics. 

In a previous study related to adjoint dynamics, 
Bcrtini et al. succeeded in deriving the large deviation 



functional S of the density profile for a special model 
of the hydrodynamic equation d t p — T>(p), which is ob- 
tained as the continuum limit for a nonequilibrium lattice 
gas |l8j. In their analysis, the equation describing the ad- 
joint dynamics dtp = T^*(p) was rigorously derived in the 
form 

2>*(P) = (x(p)V^) -A(p) (41) 

for the case in which T>(p) is given by 

V{p)= 1 -v(x{pW 5 £)+A{p), (42) 

where x(p) is the current noise intensity. They called the 
relation V(p) + V*{p) = V (x(p)V{SS/Sp)) the "fluctua- 
tion dissipation relationship" for NESSs far from equilib- 
rium, because it reduces to one expression among the lin- 
ear response relations for states near equilibrium, where 
the relation T>(p) = T>*(p) holds due to detailed balance. 
We remark that and (|4*2"j) are similar to and ((4Tj|) 
in the model we study. 

Conclusion: In conclusion, wc have found that the de- 
composition condition (|12fl leads to the relation T = fi^ , 
which was obtained in our previous study |9j. The con- 
dition (|12f> . which is related to time reversal symmetry, 
yields the new inequality l|29|l and leads to an interesting 
question regarding the decomposition of work given in 
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